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ABSTRACT 

The  purpose  of  this  research  was  to  explore  the  feasibility  of  using  spectral  hole-burning  (SHE)  memories  for 
caching  and  processing  large  amounts  of  high-speed  imagery  data  that  are  expected  in  future  space-based 
surveillance  systems.  The  research  was  to  consist  of  three  parts.  First,  the  power  of  spectral  holebuming  for  image- 
based  storage  and  processing  was  to  be  illustrated  by  performing  several  demonstration  experiments,  using  Raman 
excited  spin  coherence  (RESC).  Second,  the  properties  of  Raman  excited  spin  coherences  were  to  be  used  to 
increase  the  operating  temperature  of  spectral  holebuming  materials.  Third,  long  term  storage  of  very  high  volume 
of  optical  data  was  to  be  demonstrated  via  frequency  domain  spectral  holebuming  in  an  organic  dye  embedded  in  a 
host  of  PMMA.  In  pursuing  application  of  RESC,  we  have  demonstrated  that  nitrogen-vacancy  color  centers  in 
diamond  is  a  viable  medium  for  SHE  memory.  Furthermore,  we  have  demonstrated,  for  the  first  time  in  a  solid,  the 
slowing  and  halting  of  light  pulses  in  a  crystal  of  PriYSO.  This  effect  has  spawned  a  wide  range  of  activities  in  the 
research  community,  and  has  a  range  of  application  including  optical  data  buffering  and  high-fidelity  quantum 
memory.  For  high  volume  storage  systems,  it  became  clear  early  on  that  the  objectives  of  this  project  could  be 
better  met  by  simple  variations  of  the  tasks  originally  envisioned.  Specifically,  we  realized  that  an  organic  dye  called 
Phenanthraquinone  doped  in  PMMA  could  achieve  the  same  storage  capacity  as  envisioned  for  the  SHE  medium, 
with  a  key  advantage  that  it  could  be  used  at  room  temperature.  As  such,  we  pursued  the  development  of  this 
technology  ardently  under  this  program,  and  developed  a  novel  device  called  the  super-parallel  holographic  optical 
correlator  that  is  likely  to  be  of  near  term  practical  usage  in  space-borne  applications  such  as  missile  recognition  and 
tracking.  We  have  also  shown  how  this  device  can  be  operated  as  a  high-speed  holographic  RAM  when  operated  in 
reverse.  We  have  also  shown  how  holograms  based  on  this  material  can  be  used  to  realize  a  high  speed  holographic 
stokesmeter,  which  in  turn  can  be  employed  for  polarization  imaging  capable  of  discriminating  between  plumes  and 
rocket  bodies  during  the  boost  phase  of  a  rocket  launch.  The  device  can  also  perform  shape  discrimination  without 
imaging. 
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V  II 


2.  REPORT  OF  INVENTIONS 

•  Demonstrated  a  suepr-parallel  holographic  optical  correlator  for  ultrafast  database  search 
that  can  be  used for  missile  recognition  and  tracking 

•  Demosntrated  a  super-parallel  holographic  RAM  for  high  speed  data  recall 

•  Demonstrated  s  holographic  stokes-meter  for  polarization  imaging  that  can  be  used for 
imaging  through  plumes,  as  well  as  for  shape  discrimination  without  imaging 

•  Demonstrated  slowing  and  halting  of  light  pulses  using  Raman  excited  spin  echo  in  Pr:  YSO, 
with  application  to  data  buffering  for  optical  communication,  as  well  as  for  high  fidelity 
quantum  memory 

•  Demonstrated  the  feasibility  of  using  nitrogen-vacancy  color  center  in  diamond  as  a  viable 
medium  for  SHE  memory 
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3.  SCIENTIFIC  PROGRESS  AND  ACCOMPLISHMENTS 

The  purpose  of  this  research  was  to  explore  the  feasibility  of  using  spectral  hole-burning  (SHE) 
memories  for  caching  and  processing  large  amounts  of  high-speed  imagery  data  that  are 
expected  in  future  space-based  surveillance  systems.  The  research  was  to  consist  of  three  parts. 
First,  the  power  of  spectral  holebuming  for  image-based  storage  and  processing  was  to  be 
illustrated  by  performing  several  demonstration  experiments,  using  Raman  excited  spin 
coherence  (RESC).  Second,  the  properties  of  Raman  excited  spin  coherences  were  to  be  used  to 
increase  the  operating  temperature  of  spectral  holebuming  materials.  Third,  long  term  storage  of 
very  high  volume  of  optical  data  was  to  be  demonstrated  via  frequency  domain  spectral 
holebuming  in  an  organic  dye  embedded  in  a  host  of  PMMA.  In  pursuing  application  of  RESC, 
we  have  demonstrated  that  nitrogen-vacancy  color  centers  in  diamond  is  a  viable  medium  for 
SHE  memory.  Furthermore,  we  have  demonstrated,  for  the  first  time  in  a  solid,  the  slowing  and 
halting  of  light  pulses  in  a  crystal  of  Pr;YSO.  This  effect  has  spawned  a  wide  range  of  activities 
in  the  research  community,  and  has  a  range  of  application  including  optical  data  buffering  and 
high-fidelity  quantum  memory.  For  high  volume  storage  systems,  it  became  clear  early  on  that 
the  objeetives  of  this  project  could  be  better  met  by  simple  variations  of  the  tasks  originally 
envisioned.  Specifically,  we  realized  that  an  organic  dye  called  Phenanthraquinone  doped  in 
PMMA  could  achieve  the  same  storage  capacity  as  envisioned  for  the  SHE  medium,  with  a  key 
advantage  that  it  could  be  used  at  room  temperature.  As  such,  we  pursued  the  development  of 
this  technology  ardently  under  this  program,  and  developed  a  novel  device  called  the  super¬ 
parallel  holographic  optical  correlator  that  is  likely  to  be  of  near  term  practical  usage  in  space- 
borne  applications  such  as  missile  recognition  and  tracking.  We  have  also  shown  how  this 
device  can  be  operated  as  a  high-speed  holographic  RAM  when  operated  in  reverse.  We  have 
also  shown  how  holograms  based  on  this  material  can  be  used  to  realize  a  high  speed 
holographic  stokesmeter,  which  in  turn  can  be  employed  for  polarization  imaging  capable  of 
discriminating  between  plumes  and  rocket  bodies  during  the  boost  phase  of  a  rocket  launch.  The 
device  can  also  perform  shape  discrimination  without  imaging. 

3.1  The  Super-Parallel  Architecture  For  Holographic  Correlation  And  A 
Holographic  RAM 

[see  attachment  1  for  further  details] 

It  has  recently  been  shown  that  a  quantum  computer  could  be  used  to  factor  large  numbers 
efficiently  using  Shor's  algorithm  (SA)',  and  to  search  a  database  of  N  elements  in  steps 
using  Grover's  algorithm  (GA)^.  While  SA  offers  an  exponential  speedup  over  the  known 
classical  algorithm  for  factoring,  it  is  yet  to  be  proven  that  there  exists  no  classical  algorithm  that 
is  as  fast  as  SA.  In  contrast,  it  has  been  claimed  that  no  classical  technique  can  perform  a  search 
as  fast  as  GA.  This  claim  is  inaccurate,  based  on  the  limiting  assumption  that  a  classical  search 
would  be  carried  out  in  a  serial  manner.  Several  groups  have  shown  how  GA  does  not  rely  on 
entanglement,  but  rather  is  rooted  in  the  parallelism  of  wave  mechanics.  Such  a  parallelism  is 
also  present  in  classical  optics.  As  such,  it  is  eminently  plausible  that  a  quadratic  speed-up  in 
searching  an  imsorted  database  may  be  implementable  using  techniques  of  classical  optics.  * 
Zubairy  et  al.  recently  proposed  a  quasi-classical  technique  that  can  perform  a  search  as  fast  as 
GA,  using  unentangled  atoms,  and  parallelism  of  optics^.  Recently,  we  have  demonstrated 
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experimentally  the  basic  features  of  a  classical  search  mechanism  that  is  as  fast  as 
Specifically,  this  technique  is  based  on  super-parallel  holography  (SPH),  which  is  a  variant  of  the 
conventional  holographic  optical  correlator  (HOC)’.  While  the  speed  of  GA  scales  the  same  way 
as  that  of  these  classical  and  quasi-classical  techniques,  there  is  one  potentially  significant 
advantage  in  using  GA.  Specifically,  a  quantum  computer  needed  to  implement  a  GA  based 
search  of  an  N-element  database  requires  0(log2N)  quantum  bits,  while  any  classical  or  quasi- 
classical  technique  requires  0(N)  bits.  However,  this  seemingly  significant  advantage  of  the  GA 
is  not  of  practical  interest,  because  a  database  searching  process  is  likely  to  require  a  storage 
system  with  a  number  of  memory  bits  0(N),  even  when  the  searching  process  itself  uses  a 
quantum  computer.  Thus,  for  any  search  technique,  including  GA,  the  net  resources  necessary 
is  0(N).  Thus,  given  a  choice  of  search  engines  each  of  which  takes  steps,  and  employs 
resources  0(N),  the  preference  must  be  based  on  engineering  considerations  such  as  reliability 
and  practicability.  In  this  context,  the  SPH-based  search  engine  mentioned  above  may  be  the 
preferred  one,  given  that  it  may  be  realized  for  practical  applications  within  a  few  years,  for  non¬ 
trivial  applications  such  as  face  recognition  or  computer  vision. 

The  SPH-based  optical  correlator  (SPHOC),  invented  under  this  program'*'^,  is  a  device 
that  performs  the  holographic  optical  correlation  (HOC)  operation  simultaneously  at  thousands 
of  locations,  with  essentially  the  same  number  of  detector  elements  as  needed  for  a  single  HOC. 
To  see  how  it  performs  a  search  in  steps,  consider  first  the  operation  of  the  HOC,  which  is 
based  on  M=  images  written  in  a  single  location,  via  angular  multiplexing,  through 
interference  with  a  plane  wave  reference  beam.  When  the  query  image  is  applied  to  this 
location,  a  plane  wave  is  produced  at  the  angle  corresponding  to  the  matching  image, 
representing  a  correlation.  A  detector  array  containing  M  elements  can  be  used  to  determine  the 
identity  of  the  image.  Consider  next  a  situation  where  M  different  spatial  locations  are  used, 
each  storing  M  images.  Suppose  also  that  the  query  image  is  first  cloned  into  M  copies^,  and  the 
correlation  is  carried  out  simultaneously  at  each  of  the  M  locations.  In  a  simple  configuration, 
one  must  still  use  M  different  detector  arrays,  each  containing  M  different  detector  elements. 
The  total  number  of  detector  elements  one  must  sort  through  is  M’=N,  so  that  the  search  still 
takes  0(N)  steps.  The  SPHOC  architecture  starts  with  this  model,  using  M-fold  cloning  of  the 
query  image.  However,  a  holographic  demultiplexing  scheme  (HDS)  is  used  at  the  output  of  the 
correlator,  in  a  manner  so  that  only  two  detector  arrays,  each  containing  M  detector  elements,  are 
necessary:  one  to  identify  the  location,  and  the  other  to  identify  the  angle,  of  the  matched  image. 
Using  a  holographic  memory  unit  (HMU)  which  could  hold  two  terabyte  (TB)  of  data®,  in  a  CD- 
ROM  sized,  1  cm-thick  substrate,  the  SPHOC  searching  speed  may  be  as  high  as  one  million 
times  faster  than  a  Digital  Signal  Processing  (DSP)  based  system*. 

As  mentioned,  we  have  already  demonstrated  the  feasibility  of  the  SPHOC  on  a  small 
scale.  However,  in  order  to  realize  the  SPHOC  on  a  non-trivial  scale,  the  HDS  requires  a  large 
array  of  microscopic,  special-purpose  optical  components  that  can  be  created  only  using 
lithographic  integrated  optics  technology.  Under  this  project,  we  have  developed  a  novel 
approach  for  fabricating  a  custom-designed,  inhomogeneous,  prism-coupled  lenslet  array  for 
the  HDS  stage  of  the  architecture.  We  have  analyzed  the  actual  beam  propagation,  and 
calculated  the  exact  design  parameters  for  the  prism-lenslet-array.  In  a  future  endeavor,  the  array 
would  be  fabricated,  and  incorporated  into  a  large  scale  version  of  the  SPHOC. 


THE  SUPER  PARALLEL  HOLOGRAPHIC  OPTICAL  CORRELATOR; 
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The  key  feature  of  the  proposed  work  is  the  SPHOC  architecture,  and  variation  thereof  in  the 
form  of  an  SPHRAM.  Here,  we  describe  briefly  the  SPHOC  architecture,  shown  schematically  in 
figure  1.  First,  rXs  number  of  images  are  stored  in  each  spot  of  the  HMU  via  two-dimensional 
angular  multiplexing,  r  being  the  number  of  horizontally  multiplexed  images,  and  s  that  of 
vertically  multiplexed  images.  The  proeess  starts  by  expanding  the  collimated  beam  from  a  20 
mW  solid-state  green  laser,  for  example,  in  order  to  match  the  size  of  the  spatial  light  modulator 
(SLM).  The  target  image  is  gathered,  for  example,  by  a  camera  or  an  SAR  (synthetic  aperture 
radar).  Alternatively,  the  target  could  represent  an  object-to-image  data  map,  which  requires  a 
much  smaller  amount  of  information  to  represent  the  image.  A  copy  of  the  image  is  sent  to  the 
control  computer  for  recording,  and  another  copy  is  sent  directly  to  the  SLM  via  a  high-speed 
data  bus.  The  beam  reflected  from  the  face  of  the  SLM  carries  this  image.  This  reflected  beam  is 
now  passed  through  the  IFBR  (Image  Flattening  Beam  Reducer),  which  is  a  combination  of  two 
lenses.  It  reduces  the  image  size  by  the  ratio  of  the  focal  lengths  of  the  two  lenses,  and  projects  a 
very  slowly  diverging  version  of  the  SLM  image  at  the  face  of  the  next  element:  the  HDMX 
(Holographic  Demultiplexer/Multiplexer).  The  HDMX  is  a  cloning  device  that  produces  nXm 
copies  (which  is  the  number  of  spatial  locations  on  the  HMU)  of  the  input  image,  in  as  many 
angular  directions.  Such  a  device  can  be  constructed  simply  by  writing  nXm  plain  wave 
holograms  over  the  entire  volume. 

The  output  beams  from  the  HDMX  are  now  passed  through  an  element  termed  the  HR: 
the  Holographic  Redirector.  This  role  of  this  device  is  to  redirect  the  output  beam  from  the 
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images  now  impinge  on  the  HMU,  one  at  each  spatial  location.  If  the  image  matches  one  in  one 
of  these  locations,  a  correlation  diffracted  beam  will  emanate  from  that  location  in  a  particular 
angle.  Here,  we  show  all  such  diffraction  beams  (in  one  dimension,  for  simplicity  of  Rawing), 
keeping  in  mind  that  only  one  of  these  beams  (shown  by  the  red/bold  line)  will  be  produced  for  a 
given  image.  The  goal  of  the  rest  of  the  architecture  is  to  identify  (a)  the  spatial  location,  and  (b) 
the  angle  the  beam  corresponds  to.  These  two  pieces  of  information  will  identify  uniquely  which 
image  we  have  matched. 

The  beams  coming  out  of  the  HMU  need  to  be  collimated  in  the  direction  normal  to  the 
HMU.  This  can  be  achieved  via  the  use  of  a  prism-coupled  lenslet  array  (PCLA).  The  PCLA 
consists  of  two  stages.  The  first  stage  is  a  simple,  homogeneous,  macroscopic  lens  array 
(HMLA),  each  of  which  has  a  diameter  as  big  as  the  corresponding  storage  location  of  the  HMU. 
The  HMLA  is  to  be  positioned  abutting  the  HMU.  The  second  stage  of  the  PCLA  is  a 
microscopic,  inhomogeneous  array  of  diffractive  optic  elements  (lADOE),  each  of  which 
consists  of  a  prism,  a  lens,  and  a  phase-shifter,  all  in  series.  The  lADOE  will  be  placed  at  the 
focal  plane  of  the  HMLA.  By  "inhomogenous"  we  mean  that  each  set  of  these  prism-lens- 
phaseshifter  combination  will  have  unique  characteristics  corresponding  to  its  location.  To 
understand  how  this  would  work,  consider  a  particular  HMLA  element.  The  different  plane 
wave  (in  practice,  Gaussian)  correlation  beams  generated  from  the  HMU  area  in  front  of  it  will 
impinge  on  the  HMLA  element  at  different  angles.  As  such,  they  will  all  focus  on  the  plane  of 
the  lADOE.  The  phase  front  of  eaeh  of  these  beams  will  have  three  factors:  a  constant  phase 
shift,  a  linearly  varying  phase  shift,  and  a  quadratic  phase  shift.  The  quadratic  phase  shift  will  be 
the  same  for  each  of  these  beams;  however,  the  constant  phase  shift  and  the  linearly-varying 
phaseshifts  will  depend  on  the  angle  at  which  these  beams  were  incident  on  the  HMLA  element. 
For  the  rest  of  the  archiotecture  to  work,  we  would  like  each  of  these  beams  to  emerge  in 
parallel,  with  a  Rayleigh  length  large  enough  to  accommodate  the  rest  of  the  elements  to  follow. 
The  individual  elements  of  the  lADOE  will  be  designed  to  compensate  for  these  varying  phase 
factors,  so  that  after  exiting  the  lADOE,  all  the  beams  appear  to  be  essentially  identical,  expect 
for  a  translational  shift.  The  lens  in  the  lADOE  will  be  used  to  collimate  each  beam  ,  the  prism 
(in  the  lADOE)  will  be  used  to  compensate  for  the  linearly  varying  phase  shift,  and  the  phase- 
shifter  (in  the  lADOE)  will  be  used  to  compensate  for  the  constant  phase  shift. 

The  beams  emerging  from  the  PCLA  are  now  split  into  two  components  by  a  beam 
splitter.  One  component  goes  through  another  lenslet  array,  thus  getting  refocused,  onto  an  nXm 
array  of  CCD  elements  (CCD A).  Identifying  the  detector  that  sees  the  highest  signal  yield  the 
information  about  which  spatial  location  the  matching-image  is  in.  In  practical  situations,  the 
location-finder  array  may  suffer  from  low  signal-to-background  contrast  (STBC).  However, 
with  a  suitably  chosen  signal  integration  time  and  a  threshold  value,  it  should  still  be  possible  to 
identify  the  location  with  any  desirable  degree  of  certainty.  Alternatively,  (or  in  addition,)  the 
STBC  can  be  enhanced  by  adding  an  array  of  MEMS  based  reflectors  (controlled  by  the  angle- 
indicator  array,  as  described  in  the  next  section),  followed  by  a  low-pass  spatial  filter. 

The  other  output  of  the  beam  splitter  goes  through  another  HR  (Holographic  Redirector). 
The  HR  simply  redirects  all  the  incident  batches  to  a  central  point,  without  focusing  the  beams 
coming  from  the  same  spot  with  respect  to  one  another,  as  shown.  These  beams  are  now  passed 
through  another  HMDX,  which  is  identical  to  the  one  used  at  the  input,  but  now  operating  in 
reverse.  However,  the  reverse  operation  has  the  potential  problem  that  additional  beam  patterns 
(weaker  than  the  one  generated  along  the  axis)  will  also  be  produced.  A  simple  aperture  can  be 
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used  to  eliminated  these  unwanted  spots.  After  the  filter,  we  use  a  beam  expander  to  match  the 
size  of  the  second  CCD  array,  which  contains  rXs  number  of  elements,  corresponding  to  the 
number  of  angles  used  in  the  two-dimensional  angular  multiplexing  during  the  writing  stage.  The 
element  of  the  CCD  that  sees  the  brightest  signal  yields  the  information  about  the  angle  of  the 
matched  image.  Data  gathered  from  the  detector  array  and  the  CCD  array,  properly  thresholded, 
can  be  sent  through  a  digital  logic  circuit  to  yield  the  absolute  identity  of  the  image  we  have 
matched. 

DEMONSTRATING  THE  CONCEPT  OF  SPHOC: 

Thick  holograms  have  the  ability  to  store  multiple  images  (multiplex)  in  the  same  volume.  The 
images  can  be  multiplexed  in  several  ways.  Figure  2  demonstrates  a  technique  known  as  angle 
multiplexing.  The  holograms  are  written  by  exposing  the  substrate  with  the  image  beam  and  a 
reference  beam  at  a  specific  angle  (Figure  2(a.)).  A  second  image  is  then  written  with  the 
reference  beam  positioned  at  a  different  angle  (Figure  2(b.)).  In  this  manner,  thousands  of  images 
can  be  stored  (Figure  2(c.)).  The  images  are  searched  by  illuminating  the  sample  with  a  query 
image  (Figure  2(d.)).  If  the  image  matches  with  a  stored  image,  a  diffracted  beam  emerges  from 
the  back  of  the  sample.  The  direction  of  this  beam  indicates  which  image  has  matched. 
Therefore,  all  of  the  images  are  searched  simultaneously. 

Figure  3  displays  images  taken  directly  from  the  spatial  light  modulator  (SLM)  as  well  as 
images  that  were  read  from  the  HMU.  The  images  are  read  from  the  HMU  by  illuminating  the 
sample  with  the  appropriately  angled  reference  beam.  A  Cohu  CCD  camera  was  used  to  capture 
the  images,  and  transferred  to  the  computer  using  a  frame-grabbing  software.  The  columns 
marked  "SLM  IMAGE"  shows  images  that  are  produced  directly  by  the  SLM,  and  captured  by 
the  CCD  camera.  The  columns  marked  "HOLOGRAPHIC  IMAGE"  show  images  that  were  read 
out  from  the  holograms,  and  then  captured  by  the  same  CCD  camera 
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Figure  2..  Angle-multiplexed  holographic  image  storage 
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Figure  3 :  Sample  images  stored  and  recalled  from  the  HMU. 


Figure  4.  Eight  images  stored  in  one  spatial  location  of  a 
HMU  unit.  Each  image  was  produced  by  illuminating  the 
hologram  with  a  plane-wave  beam  at  the  appropriate 
angle. 


As  can  be  seen,  reasonably  high-quality  images  were  successfully  stored  in  the  PDA 
(photopolymer  with  diffiision  amplification)  material  used  for  this  project.  The  PDA  material 
was  invented  originally  at  the  Vavilov  State  Optical  Institute  (VSOI)in  St.  Petersburg,  Russia, 
and  its  optical  quality  was  enhanced  further  through  collaborative  work  carried  out  between 
VSOI,  the  group  of  the  PI,  and  the  Laser  Photonics  Technology  (LPT)  of  Amherst,  NY.  We  can 
now  produce  this  material,  in  virtually  any  custom  size  and  shape,  in  our  laboratory  at  NWU. 
Furthermore,  LPT  markets  a  variant  of  this  material  under  the  trade  name  of  Memplex™.  These 
holograms  represent  the  best  quality  holographic  images  we  have  recorded  to  date.  In  order  to 
achieve  these  results,  we  had  to  make  several  improvements  to  our  setup,  including  robust 
vibration  isolation,  and  switching  from  a  rather  unstable  Ar^  laser  to  the  super-stable  Verdi  laser. 
Figure  4  shows  another  set  of  images  stored  in  a  single  location  of  an  HMU  via  angular 
multiplexing. 

Using  the  same  exposure  and  curing  parameters  which  had  previously  resulted  in  high 
quality  stored  holographic  images  in  the  HMU,  we  created  an  HMU  with  three  spatial  locations 
by  translating  the  sample  in  one  dimension.  The  memory  cells  were  placed  0.25”  apart  in  order 
to  avoid  cross-talk  between  the  cells.  The  contents  of  each  cell  were  verified  by  illuminating  the 
sample  with  reference  beams  and  observing  the  holographically  re-generated  images.  This  was  a 
key  step  in  verifying  the  scalable  nature  of  the  SPHOC  architecture.  Following  our  successful 
writing  of  a  HMU  unit,  we  next  attempted  to  correlate  the  stored  images  with  query  images. 


In  order  to  facilitate  good  spatial  alignment  of  the  HMU  with  the  query  beam,  new  stored 
images  were  composed.  The  images  consisted  of  a  three  bladed  “fan”  that  is  rotated  by  15 
degrees  with  respect  to  the  preceding  image.  In  this  way  there  should  be  no  spatial  overlap 
between  the  test  images.  An  example  image  is  shown  in  Figure  5.  Using  this  new  design  another 
HMU  unit  was  written  using  the  same  exposure  and  cure  parameters  used  to  produce  the 
previous  sample.  We  then  tested  the  sample  to  insure  that  the  holographic  images  had  been 
stored.  We  next  precisely  aligned  the  sample  with  the  original  SLM  image  while  observing  the 
diffracted  rays  that  emerged  from  the  sample.  The  correlation  occurs  as  follows:  when  the  image 
projected  on  the  sample  matches  one  of  the  stored  images  a  plane  wave  emerges  in  a  specific 
direction.  We  collected  all  of  the  plane  or  diffracted  waves  with  a  large  lens  which  Fourier 
transforms  the  unique  angles  into  unique  spatial  locations.  In  this  way  the  matching  image  can  be 
determined  by  observing  the  diffraction  spot  intensity.  Figure  9  summarizes  the  results  of  this 
correlation  study.  The  images  in  Figure  6  represent  the  query  image  projected  onto  the  HMU. 
Column  II  displays  experimental  images  of  the  diffraction  spots,  and  Column  III.  displays 
quantitative  intensity  measurements  taken  along  the  entire  row  of  spots.  As  can  be  seen  from  the 
data,  the  correlation  shows  strong  peaks  with  signal  to  noise  ratios  on  the  order  of  150.  This  data 
validates  the  ability  of  the  thick  holographic  correlator  to  discriminate  images  clearly. 


Figure  5:  An  example  of  the  type  of 
images  used  for  demonstrating 
efficient  correlation.  This  image  was 
read  from  the  HMU. 


Figure  6.  Demonstration  of  HMU  image  correlation 


In  order  to  demonstrate  parallel  image  correlation,  the  query  image  was  directed  through 
a  1x3  HMDX  or  “splitter”.  This  holographic  optical  element  produces  three  copies  of  the 
incident  image.  The  original  SLM  image  as  well  as  the  resulting  three  projected  images  are 
displayed  in  Figure  7.  The  HMDX  allows  the  query  image  to  be  projected  simultaneously  on 
every  memory  cell.  In  order  to  test  each  optical  element,  the  split  images  were  then  projected 
onto  the  three-cell  HMU  with  the  use  of  mirrors.  This  provided  us  with  the  flexibility  to  direct 
each  beam  individually,  making  query-beam/HMU  alignment  slightly  easier.  The  diffracted 
beams  generated  by  the  HMU  were  then  collected  by  a  lenslet  array  (shown  in  fig  8)  and 
projected  onto  a  CCD  camera  for  quantitative  intensity  data.  A  schematic  of  this  set-up  is 
included  in  Figure  9. 
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Figure  7.  Image  on  the  left  is  the  SLM  generated  test  image.  On  the  formes  by  a  commerically  available  8X8  lenslet 
ri^t  is  an  experimental  image  of  the  output  of  the  1x3  HMDX,  array,  housed  in  a  circular  substrate. 
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Figure  9.  Schematic  of  correlation  set-up 


Figure  10.  Correlation  data  from  three  multiplexed  query  imagps 


The  collected  data  is  displayed  in  Figure  10.  Each  column  represents  output  from  each 
HMU  memory  cell.  The  images  in  the  far-left  column  represent  the  query  images  presented  to 
the  HMU  unit.  The  same  images  were  stored  in  the  HMU.  The  images  in  columns  I,  II,  and  HI 
display  the  CCD  images  of  the  diffracted  beams  that  emerged  from  the  HMU  when  each  of  the 
query  beams  (I,  II,  and  III)  were  directed  onto  the  HMU.  The  lateral  position  of  the  bright  spots 
(a,  b,  c. . .)  indicates  which  of  the  stored  images  has  matched  with  the  query  image.  As  the  data  in 
Figure  10  shows,  the  HMU  is  providing  a  clear  indication  of  which  of  the  stored  image  the 
queried  image  represents.  The  signal  to  noise  ratio  of  this  data  is  quite  good,  given  the  low 
resolution  of  om  SLM. 

The  next  step  would  be  to  extend  this  architecture  to  two  dimensions.  To  this  end,  we 
have  constructed  and  demonstrated  a  1X9  HMDX,  with  the  nine  output  images  arranged  in  a 
symmetric  3X3  configuration.  The  function  of  this  HMDX  is  illustrated  in  figure  11,  with  one  of 
the  fan  images  as  a  test  input.  As  can  be  seen,  the  images  are  split  as  expected,  without  any 
significant  distortions. 
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Figure  11.  Demonstration  of  holographic  multiplexer/demultiplexer 
(HMDX);  (a)  schematic  of  the  3X3  HMDX,  (b)  SLM  image  that  was 
input  to  the  HM  DX,  and  (c)  CCD  capture  of  the  nine  resulting  images 


SUPER-PARALLEL  HOLGRAPHIC  RANDOM  ACCESS  MEMORY: 

For  a  practical  image  recognition  and  tracking  system  (IRTS),  it  is  useful  to  have  a  high-speed 
random  access  memory  (RAM)  to  complement  the  SPHOC.  However,  the  speed  of  a 
conventional  RAM,  which  accesses  data  serially,  is  too  slow  for  this  application.  Here,  we  show 
how  a  super-fast  RAM  can  be  realized  by  simply  operating  the  SPHOC  in  reverse.  Called  the 
super-parallel  holographic  random  access  memory  (SPHRAM),  this  device  and  the  SPHOC 
can  be  operated  alternately  using  essentially  the  same  hardware,  and  potentially  the  same  data 
base. 


Figure  12.  Schematic  illustration  of  the  Super-Parallel  Holographic  RAM 
(SPHRAM)  architecture 
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The  basic  idea  of  the  SPHRAM  is  illustrated  in  Figure  12.  As  stated  above,  the  functioning  of 
this  device  is  similar  to  an  SPHOC  operated  in  reverse.  Briefly,  the  holographic  memory  unit 
(HMU)  is  recorded  with  the  database  of  interest,  using  multiple  spatial  locations,  each  of  which 
contains  a  set  of  images  that  are  angularly  multiplexing  in  two  dimensions.  When  operating  the 
RAM,  the  user  enters  as  coordinates  the  spatial  position  and  the  angle  of  storage  of  the  image  of 
interest.  The  position  coordinate  is  used  to  open  the  corresponding  element  of  a  shutter 
array.  The  number  of  elements  in  the  shutter  array  is  the  same  as  the  number  of  spatial  locations 
in  the  HMU,  and  the  locations  are  matched.  Shutter  arrays  of  this  type  are  available 
commercially,  with  individually  addressable  ferroelectric  liquid  crystals,  for  example,  with  a 
very  fast  switching  time.  The  extinction  ratio  of  these  shutters  can  be  as  high  as  10^;  cascading  a 
pair  of  shutters  can  be  used  to  achieve  an  extinction  ratio  of  lO’*’. 

The  angular  coordinate  for  the  image  is  sent  to  a  beam  deflector  (e.g.,  using  a  pair  of 
acousto-optic  or  electro-optic  deflectors),  which  orients  the  read  beam  at  the  desired  angle, 
which  in  turn  is  translated  to  a  specific  position  by  the  redirector.  (Alternatively,  one  could 
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Figure  13:  Examples  of  data  pages  retrieved  through 
a  simplified  version  of  the  SPHRAM  architecture, 
using  a  pair  of  galvo-mounted  mirrors  for  deflections 

simply  use  an  SLM,  but  it  would  reduce  the  amount  of  light  available  for  read-out).  The 
combination  of  the  reducing  telescope,  the  multiplexer,  the  redirector,  and  the  lenslet  array  (in 
practice,  the  PCLA  described  above)  produces  a  copy  of  the  read  beam  simultaneously  at  each  of 
the  locations;  this  is  the  key  feature  of  the  SPHRAM.  As  such,  the  image  stored  at  this  angle 
will  be  recalled  from  each  location.  The  shutter  array  will  block  all  but  one  of  these  images,  and 
the  redirector  and  the  demultiplexer  will  send  the  wanted  image  to  the  CCD  camera. 


We  have  recently  demonstrated  the  feasibility  of  one  part  of  the  concept  here,  by  using  a 
pair  of  galvo-mounted  mirrors  for  deflections.  The  geometry  used  by  us  for  the  demonstration  is 
shown  in  figure  13.  For  the  database,  we  used  an  existing  HMU  which  has  9  locations,  in  a  3X3 
arrangement,  each  location  containing  8  images  multiplexed  angularly.  The  redirectors  and  the 
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splitter  (multiplexer)  used  are  the  same  as  the  ones  created  before  for  demonstrating  the  SPHOC 
architecture.  The  bottom  of  the  figure  shows  a  typical  set  of  data  retrieved  using  this  setup. 


Figure  14:  Schematic  illustration  of  an  architecture  that 
can  be  used  as  either  an  SPHOC  or  an  SPHRAM 

Similar  data  were  also  retrieved  from  the  other  locations  (not  shown). 

Note  that  a  simple  modification  of  this  architecture  will  allow  one  to  operate  the  same 
system  as  the  SPHOC  as  well.  This  is  illustrated  in  figure  14.  During  the  operation  as  an 
SPHOC,  all  the  elements  of  the  shutter  array  are  open,  shutter  SI  is  closed,  and  shutter  S2  is 
open.  Beams  reflected  by  the  beam  splitter  BSl  is  used  for  location  detection,  and  the  beams 
reflected  by  the  beam  splitter  BS2  is  used  for  angle  detection.  During  the  operation  as  an 
SPHRAM,  SI  is  opened,  S2  is  closed,  and  the  elements  of  the  shutter  array  is  opened  selectively 
according  to  the  address  of  the  data  page  to  be  retrieved. 

ARCHITECTURAL  CHALLENGES  ASSOCIATED  WITH  THE  SPHOC  AND  THE  SPHRAM,  AND  THE 
PROPOSED  SOLUTION:  BEAM  PROPAGATION  ANALYSIS  OF  THE  SPHOC  GEOMETRY 

One  of  the  main  problems  associated  with  the  SPHOC/SPHRAM  architecture  is  that  of  beam 
collimation  after  the  HMU.  The  diffraction  from  the  holograms  recorded  at  large  angles 
propagate  at  an  angle  with  respect  to  the  optic  axis  of  the  system.  These,  when  collimated  using 
one  of  the  lenses  of  the  HMLA  (homogeneous,  macroscopic  lenslet  array),  do  not  propagate 
parallel  to  the  optic  axis,  thus  creating  a  problem  in  the  detector  plane.  As  discussed  above  (see 
the  two  paragraphs  following  figure  1),  this  problem  can  be  solved  by  adding  the  PCLA,  to  be 
placed  at  the  focal  plane  of  the  HMLA.  Here,  we  discuss  these  issues  explicitly,  using  a  model 
where  the  correlation  beam  produced  is  a  plane  wave.  For  a  real  system,  this  analysis  has  to  be 
augmented,  using  FDTD  (finite  division  time  domain)  techniques,  by  the  facts  that  (a)  the 
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correlation  spots  would  be  Gaussian  in  nature,  and  (b)  the  exact  shape  of  the  correlation  spots 
would  depend  on  the  thick-holographic  correlation  process.  The  augmented  analysis  would  be 
one  of  the  tasks  carried  out  under  this  project.  Following  the  simplified  beam  propagation 
analysis  of  the  system,  we  present  the  lithographic  design  of  the  prism-coupled  lenslet  array 
(PCLA)  that  we  propose  to  make  at  Northwestern  University  under  this  project  proposed  here. 


PROPAGATION  ANALYSIS  FOR  THE  BEAM  PATH  AFTER  THE  HMU 

In  order  to  obtain  a  successful  detection  of  a  match  of  any  given  HMU  image,  the  SPHOC 
architecture  must  be  able  to  collimate  and  direct  every  possible  correlation  beam  towards  the 
detector  array.  For  a  large  capacity  HMU,  there  will  be  many  angles  and  thus  the  possible 
correlation  beams  will  overlap  to  a  large  extent  (Fig  15).  In  order  to  resolve  each  beam,  a  lens 
(i.e.,  an  element  of  the  HMLA)  directly  following  the  HMU  can  be  used  to  focus  each  angled 
beam  to  a  different  location  in  the  focal  plane.  However,  the  beams  will  still  be  traveling  at 
angles  to  the  optic  axis,  so  redirection  and 


Figure  15a:  Overlapping  Correlation  Beams  Figure  15b:  Lens  to  separate  beams 


We  propose  to  use  the  prism-coupled  lenslet  array  (PCLA)  described  above  (see  the  two 
paragraphs  following  figure  1)  a  short  distance  after  the  focal  plane  of  this  initial  lens.  This  will 
redirect  and  collimate  the  beams  as  well  as  de-magnify  the  beams  by  the  ratio  of  the  focal 
lengths  of  the  lens  chosen.  Here  we  use  a  simple  method  to  determine  the  proper  choice  of 
parameters  to  ensure  collimation  and  redirection  parallel  to  the  optic  axis  of  the  system. 

Consider  the  case  of  a  single  correlation  beam  traveling  at  an  angled  with  respect  to  the 
optic  axis  z  .  This  beam  will  enter  the  initial  lens  and  will  be  offset  by  a  certain  amount  in  the 
focal  plane  (Fig  16a). 


Figure  16a:  Single  Correlation  Beam.  Although  not  obvious  in  this  drawing  ,  after  focusing,  the  beam  will  travel  at  an  angle  with  respect 
to  the  optic  axis,  due  to  the  linearly  vaiying  phase  shift  imparted,  the  amount  of  this  being  proportional  to  the  distance  of  the  incidence 
spot  away  from  the  center  of  the  lens. 

Although  this  picture  is  accurate,  a  more  useful  way  to  view  this  is  to  consider  a  shifted  optic 
axis  z' (Fig  16b),  which  shows  explicitly  that  a  beam  incident  on  a  lens  at  a  position  Xq  away 
from  the  center  sees  a  phase  shift  corresponding  to  that  of  a  lens  center  around  the  point  of 
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incidence,  and  a  prism,  the  gradient  of  which  is  proportional  to  Xo.  Specifically,  the  following 
phase  transformation  will  be  applied  to  the  beam: 
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Thus  this  beam  is  affected  by  a  quadratic  phase,  which  focuses  the  beam,  by  a  linear  phase, 
which  redirects  the  beam,  and  by  a  constant  phase,  which  we  make  note  of  but  disregard  for  now 
(this  factor  can  also  be  compensated  for  in  the  PCLA,  as  discussed  before). 

We  can  express  a  plane  wave  traveling  at  angled  hy  U  =  exp(i2;rUgx)  where 
Vq  =  sin  OjX .  Therefore,  the  field  directly  after  the  lens-prism  combination  is: 
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We  use  the  Fresnel  integral  under  the  paraxial  approximation  to  determine  the  field  at  z=f  and 
z=f+L: 
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where  x' is  the  coordinate  in  the  z=f  plane,  x  is  the  coordinate  in  the  z=0  (initial  lens)  plane,  and 
f'  =  x'IXf.  Again  we  have  noted  but  disregarded  the  constant  phase  terms.  Evaluating  this 
integral  gives  us  the  expression  for  the  field  at  z=f: 
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We  can  now  use  the  Fresnel  integral  a  second  time  to  propagate  this  field  to  z=f+L: 


t^z=/-n=exp 


(.k^ 

r»  f .  /r  ,2^0 

/ 

1 

1 

J 

U/" 

r  2;r/JJ 

exp 

V 

jlTif”  x'\dx' 

2/  J 


In  this  expression  x"  is  the  coordinate  in  the  z=f+L  plane.  Evaluating  this  integral  gives  us  our 
final  expression  for  the  field  in  the  z=f+L  plane.  We  note  that  it  has  a  quadratic  phase  factor  as 
expected  by  the  focusing  action  of  the  lens,  and  that  the  field  is  traveling  tilted  with  respect  to  the 
optic  axis. 
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This  yields  the  information  on  the  focal  length  of  the  lens  (in  the  PCLA)  we  will  need  to  cancel 
the  quadratic  phase  term,  as  well  as  the  angle  of  the  prism  (also  in  the  PCLA,  juxtaposed  to  the 
lens)  needed  to  cancel  the  linear  phase  term,  thus  producing  a  plane  wave  traveling  parallel  to  the 
optic  axis.  By  keeping  track  of  all  the  constant  phase  terms  added  by  the  prisms,  lenses,  and 
free-space  propagation,  we  can  also  add  a  phase  plate  to  the  PCLA.  The  final  array  will  consist 
of  a  quadratic  phase  correction,  a  linear  phase  correction,  and  a  phase-compensator  plate, 
yielding  at  its  output  collimated  plane  waves  traveling  along  the  optic  axis,  in  phase,  as  required 
by  the  SPHOC/SPHRAM  architecture.. 

Assuming  a  3mm  beam  diameter  entering  the  first  lens  (f=lcm)  and  a  desired  diameter  of 
300  fxm ,  the  second  lens  will  have  a  focal  length  of  2mm  and  the  output  beam  will  have  a 
Rayleigh  range  of  approximately  53  cm  at  532nm.  Such  a  long  Rayleigh  range  is  more  than 
adequate  for  the  architecture  of  the  SPHOC  and  the  SPHRAM. 

LITHOGRAPHIC  DESIGN  AND  CONSTRUCTION  OF  THE  PRISM-COUPLED  LENSLET  ARRAY  AND 
OTHER  CHALLENGING  COMPONENTS: 

The  beam  path  from  the  HMU  to  the  prism-coupled  lenslet  array  is  shown  in  fig.  17  below  (the 
unfilled  dots  in  the  focal  plane  are  meant  to  all  lie  in  plane  perpendicular  to  the  propagation  axis, 
and  each  dot  represents  an  element  of  the  PCLA).  As  mentioned  earliar,  the  HMU  generates  rXs 
diffracted  beams  from  each  of  its  nXm  locations.  For  simplicity,  we  show  an  off-axis 
propagation  here.  All  the  diffracted  beams  are  focused  in  the  focal  plane  of  the  lens  (which  is  an 
element  of  the  HMLA).  The  beams  then  propagate  further  as  per  their  corresponding  divergence. 
Each  of  these  beams  is  then  collected  by  an  element  of  the  prism-coupled  lenslet  array  (PCLA), 
which  collimates  the  beam  along  the  optic  axis  of  the  system.  As  illustrated  in  detail  above,  this 
PCLA  is  to  be  custom-designed  so  that:  (a)  each  element  of  the  array  has  a  prism  with  a  imique 
deflection  angle  corresponding  to  the  linearly  varying  phase  shift  experienced  by  the  input 
beam,  (b)  the  lenses  in  the  array  have  the  same  focal  length,  and  they  collimate  the  incident 
beams  along  the  optic  axis,  and  (c)  each  element  provides  a  unique  and  constant  phase  shift 
factor. 
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Figure  17:  Beam  Propagation  from  HMU  to  the  Proposed  Prism-Coupled  Lenslet  Array 


Specifically,  one  can  fabricate  the  focusing  lens 
along  with  the  prism-coupled  lenslet  array  using  a 
silicon  bench.  This  would  enable  an  easy  pre¬ 
fabricated  alignment  of  each  individual  element  in 
the  array.  The  following  figure  (figure  18)  depicts 
one  such  assembly.  The  basic  construction  of  prism- 
coupled  lenslet  array  is  based  on  very  strong  graded- 
refractive-index  (GRIN)  structure  to  achieve  focusing 
along  the  optic  axis  of  the  SPHOC  system,  which 
will  be  vertical  to  the  Silicon  bench  on  which  we  will 
fabricate  the  whole  assembly.  Phase  correction  in  this 
direction  can  be  achieved  by  refractive  index  profile 
as  shown  in  Fig  19a.  In  the  other  direction  (horizontal 
to  Si  bench),  focusing  and  phase  correction  are 
achieved  by  physical  curved  surfaces  (lens  plus 
prism)  as  shown  in  Fig.  19(b). 


Figure  18:  Assembly  showing  the  focusing  lens  along 
with  the  prism-coupled  lenslet  array 
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Vertical  Direction 

Focusing  &  Phase  Correcting  by 
parabolic  profile  plus  linear  profile 


Lateral  Direction 

Phase  Correction  Part  Focusing  Part 


(a)  (b) 

Fig  19  Precision  lens  outline:  (a)  vertical  focusing  and  phase  correction  by  multilayer 
structure,  (b)  Lateral  focusing  and  phase  correction  by  physical  surface  profile 


Suppose  refractive  index  at  the  center  of  the  GRIN  is  rio  and  the  refractive  index  on  the 
side  is  tir,  the  effective  refractive  index  difference  is  defined  by  c/«  =  -y/«o  -n\  . 
According  to  GRIN  lens  formula,  the  numerical  aperture  NA  is  about  the  same  value  as 
dn.  The  spot  size  a  (FWHM)  is  d.SXIdn.  One  limitation  of  the  conventional  GRIN  lens 
is  that  high  NA  cannot  be  achieved,  as  the  graded  refractive  index  profile  is  achieved 
via  material  diffusion  technique.  Such  technique  cannot  achieve  large  difference  in 
refractive  indices  between  the  lens  center  and  side.  Another  limitation  is  the  refractive 
index  profile  cannot  be  arbitrary.  The  superlens  is  capable  of  correcting  for  beam’s 
spatial  aberration  and  astigmatism.  It  is  a  powerful  precision  lens  that  can  be  made 
directly  on  silicon  bench  in  large  number  with  each  lens  capable  of  performing 
individual  complex  beam  profile  transformation. 


GRIN  lens  fabrication  methods  makes  use  of  By  inter-spacing  multiple-layer  thin  films 
of  two  materials  with  drastic  difference  in  refractive  index  to  achieve  a  highly 
controllable  refractive  index  variation,  one  high  refractive  index  material  with  another 
low  refractive  index  material,  we  can  get  effective  refractive  index  of  any  value  in 
between.  The  principle  of  this  technique  is  shown  in  fig  20.  When  mixing  thin  layer  of 
two  different  refractive  index  materials,  the  effective  refractive  index  of  the  mixed 
section  can  be  calculated  as: 
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where  ni  is  the  higher  refractive  index,  «2  is  the  lower  refractive  index,  Li  is  the 
thickness  of  the  higher  refi'active  index  material  and  L2  is  the  thickness  of  the  lower 
refractive  index  material. 
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(a) 


Low  index 


Fig.  20  (a)  Two-material  layered  structure  with  parabolic  graded  refractive  index  profile;  (b)  Mix  two 
materials  to  get  the  varying  effective  index  to  achieve  desired  index. 


Using  Si02/Ti02,  a  refractive  index  variation  of  1.5  to  2.4  can  be  achieved  (Fig.  21(a)). 
In  Fig.  21(b)  below,  one  can  show  the  strong  focusing  capability  of  this  dual  material 
structure  by  Finite  Difference  Time  Domain  (FDTD)  simulation. 


Fig.  21  (a)  3D  Lens  structure,  (b)  FDTD  simulation  of  precision  lens  in  vertical  direction,  with  output  to 
150nm  channel  waveguide 

The  dual  material  structure  is  able  to  provide  any  refractive  index  profile,  such  as 
parabolic  graded  profile  plus  linear  graded  profile.  The  parabolic  profile  will  provide 
focusing  or  collimating  and  the  linear  profile  will  provide  phase  correction  as  shown  in 
Fig.  19(a). 

In  lateral  direction,  focusing  and  phase  correction  are  achieved  by  physical  curved 
surfaces,  which  can  be  any  shape  because  it  is  only  defined  by  photolithographic  mask. 
As  shown  in  Fig.  22,  the  lens  part  is  providing  focusing  and  the  prism  part  provides 
phase  correction.  To  produce  the  parallel  beam  OCM,  different  phase  correction  (or 
angular  transformation)  as  shown  in  (a),  (b),  (c)  are  needed  at  different  FP  focused 
points  (see  Fig.  22(d)).  Using  superlens,  these  slightly  different  lenses  can  be  fabricated 
in  a  single  run  of  fabrication. 


Lateral  Direction 


Lateral  Direction 


Lateral  Direction 
Phase  Correction  Part 


Overview 


Fig.  22.  Lateral  focusing  plus  different  phase  correction  (a,  b,  c,  d). 
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(b)  E-beam  pattern  with  the  dose 
distribution  writing 

(c)  E-beam  resist  development 

Figure  23.  Procedure  of  the  electron-beam  writing  fabrication  of  Fresnel  lens 


The  Fresnel  lens  can  be  fabricated  using  e-beam  lithography.  The  electron-beam 
writing  procedure  of  fabrication  Fresenl  lens  is  shown  in  Fig.  23.  First,  an  electron-beam 
resist  is  spin-coated  on  a  glass  substrate  covered  with  a  transparent  conductive  layer  of 
ITO,  which  is  necessary  to  avoid  electrical  charge  buildup  during  the  writing.  In  order  to 
increase  the  diffraction  efficiency  of  Fresnel  lens,  a  blazed  Fresnel  grating  is  required  in 
the  process.  The  initial  resist  layer  should  be  slightly  thicker  than  the  final  thickness  after 
having  been  developed.  We  can  use  PMMA  (positive  type)  or  CMS  (negative  type)  e- 
beam  resist,  where  CMS  stands  for  chloromethylated  polystyrene.  The  final  thickness  is 
determined  by  working  wavelength  and  E-beam  resist  refractive  index.  Our  system 
working  wavelength  is  ^=0.532|^m  and  E-beam  resist  refractive  index  n=1.49  (PMMA). 
The  final  thickness  is  estimated  to  be  about  T  =  A/(n-l)  =  1.08//w.  Second,  the  blazed 
profile  can  be  formed  by  a  suitably  chosen  electron-dose  distribution  because  the  etching 
rate  of  the  resist  will  depend  on  the  electron  dose.  The  electron-dose  versus  etch-depth 
relationship  is  obtained  from  a  cross-sectional  scanning-electron-microscope  photograph 
of  test  sample  that  has  been  drawn  by  linearly  changing  the  dose.  After  getting  the 
electron-dose  versus  etch-depth  relationship,  one  can  scan  the  sample  at  different 
locations  of  the  Fresnel  for  a  different  predetermined  time  until  the  specified  dose  is 
given  to  obtain  the  blazed  grating  shape  shown  in  Fig.  23c.  Third,  after  exposure  the 
sample  is  carefully  developed  in  1 : 1  solution  of  isopropyl  alcohol  and  methyl  isobutyl 
ketone.  The  development  of  the  resist  would  achieve  the  desired  relief  strucmre. 
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3.2  Theory  Of  Holographic  Beam  Combining  With  Applications  To  Super-Dense 
Wavelength  Division  Multiplexing  And  Holographic  Correlators 

We  have  developed  a  mathematical  model  of  coherent  and  incoherent  beam  combination 
in  a  thick  hologram.  We  also  derive  the  formulae  relating  the  read  and  write  angles  to  the 
read  and  write  wavelengths  for  the  combiner.  Furthermore,  we  present  a  new  technique 
for  determining  the  M#,  and  establish  that  the  M#  required  for  a  coherent  combiner  is 
substantially  less  than  that  needed  for  an  incoherent  one.  These  models  are  of  direct 
importance  to  developing  key  components  of  the  SPHOC  architecture  described  above, 
as  well  as  for  potential  applications  to  super-dense  WDM  systems  based  on  holographic 
beam  combiners. 

[SEE  ATTACHMENT 2  FOR  DETAILS] 
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3.3  Demonstration  Of  Nitrogen  Vacancy  Color  Center  In  Diamond  As  A  Viable 
Medium  For  SHE  Memory 

We  have  observed  Raman-excited  spin  coherences  in  nitrogen- vacancy  diamond  color 
centers  by  means  of  nondegenerate  four-wave  mixing  and  electromagnetically  induced 
transparency.  The  maximal  absorption  suppression  was  found  to  be  17%,  which 
corresponds  to  70%  of  what  is  possible  given  the  random  geometric  orientation  of  the 
color  center  in  diamond.  This  experiment  demonstrate  the  feasibility  of  high-capacity, 
high-speed  data  storage  via  Raman  type  SHE  process  in  this  material. 


[SEE  A  TTA  CHMENT  3  FOR  DETAILS] 


3.4  Demonstration  Of  Slowing  and  Stopping  of  a  Light  Pulse  in  a  Crystal  of 
Pr:YSO  with  Applications  to  Data  Buffering  and  Quantum  Memory 

We  have  demonstrated  ultraslow  group  velocities  of  light  in  an  optically  dense  crystal  of 
Pr:YSO.  Light  speeds  as  slow  as  45  m/s  were  observed,  corresponding  to  a  group  delay 
of  66  ms.  Deceleration  and  “stopping”  or  trapping  of  the  light  pulse  was  also  observed. 
These  reductions  of  the  group  velocity  are  accomplished  by  using  a  sharp  spectral  feature 
in  absorption  and  dispersion  that  is  produced  by  resonance  Raman  excitation  of  a  ground- 
state  spin  coherence.  This  process  has  potentially  significant  applications  in  buffering  of 
high-bit-rate  data  streams  in  optical  communication,  as  well  as  in  realizing  a  high-fidelity 
quantum  memory  system. 

[SEE  A TTACHMENT  4  FOR  DETAILS] 
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3.5  A  Holographic  Stokes-meter  for  Polarization  Imaging 

Polarimetric  imaging*'^  takes  advantage  of  the  fact  that  a  given  object  emits  and  scatters 
light  in  a  unique  way  depending  on  its  polarimetric  signature.  Identifying  the  polarimetric 
signature  is  equivalent  to  identifying  the  scattered  Stokes  vector'^’^.  An  active  polarimetric 
sensor  is  used  in  applications  such  as  target  recognition,  vegetation  mapping,  pollution 
monitoring,  geological  surveys,  medical  diagnostics  and  so  on.^‘^ 

Current  polarization  imaging  systems  include  mechanical  quarter-wave 
plate/linear  polarizer  combinations photo-detectors  with  polarization  filtering  gratings 
etched  onto  the  pixels,  and  liquid  crystal  variable  retarders.  The  speed  of  the  mechanical 
sensor  is  limited  because  each  Stokes  parameter  is  determined  sequentially  and  the  wave 
plate/polarizer  combination  must  be  re-oriented  precisely  between  parameters.  The 
etched  photo-detector  systems  cannot  resolve  the  complete  Stokes  vector  at  this  time. 
The  liquid  crystal  retardation  is  very  similar  to  the  mechanical  sensor  but  with  a  liquid 
crystal  display  replacing  the  wave  plates  and  polarizers.  This  method  is  still  sequential 
and  is  restricted  by  the  time  it  takes  the  display  to  re-orient  itself,  typically  on  the  order  of 
100ms  per  scan.  This  limits  the  throughput  to  ~  10  Hz.  To  examine  the  advantages  of 
the  proposed  architecture  and  to  quantify  its  speed,  we  examine  each  component  (see 
Figure  1).  A  typical  thick  hologram  (~lmm)  has  a  channel  bandwidth  on  the  order  of  1 
nm  (and  an  angular  bandwidth  of  ~1  mrad),  corresponding  to  an  optical  response  time  of 
~10  ps.  The  signal  manipulation  can  be  accomplished  with  pre-calibrated  field- 
programmable  gate  arrays  (FPGA)  or  programmable  logic  arrays  (PLA)  and  does  not 
require  real-time  processing.  These  devices  perform  at  roughly  the  speed  of  the  logic 
gates,  typically  on  the  order  of  1  ns.  The  detector  array  response  time  is  determined  in 

part  by  the  desired  signal  to  noise  ratio  (SNR).  The  SNR  is  proportional  to^jrjlr  ,  where 
7  is  the  quantum  efficiency,  I  is  the  intensity,  and  r  is  the  average  time.  For  a  typical 
SNR  of  1 0,  the  response  time  is  then  given  by  the  detector-specific  parameters  and  the 
amount  of  light  reflected  by  the  target.  It  should  be  noted  that  the  constraints  imposed  by 
a  given  SNR  and  the  detector  array  are  common  to  all  of  the  methods  discussed  here. 
The  advantage  of  our  design  is  that  it  can  determine  the  complete  Stokes  vector  in 
parallel,  so  that  the  detector  and  not  the  sensor  is  the  limiting  factor  in  determining  the 
speed  of  the  device. 

The  architecture  is  shown  in  Figure  1 .  The  incoming  image  is  split  into  two  copies 
using  a  beam  splitter.  The  first  beam  is  diffracted  into  two  beams  using  two  multiplexed 
holographic  gratings.  The  second  beam  passes  through  a  quarter-wave  plate  before 
diffracting  from  a  similar  set  of  multiplexed  holographic  gratings.  The  diffracted  beams 
are  projected  onto  four  CCD  arrays.  Their  intensities  are  summed  with  pre-determined 
weights  to  compute  each  component  of  the  Stokes  vector.  These  weighting  factors  are 
determined  by  using  the  parallel  and  perpendicular  polarization  components  of  diffraction 
to  formulate  the  Mueller  matrix  that  describes  the  transformation  of  the  initial  Stokes 
parameters  by  grating  diffraction.  This  design  takes  advantage  of  the  fact  that  a 
hologram  is  sensitive  to  the  polarization  of  the  incident  light,  and  the  weighting  factors 
can  be  determined  analytically  by  using  a  Mueller  matrix  analysis  of  the  architecture. 

For  an  arbitrary  image  pattern,  the  diffraction  efficiencies  depend  also  on  the 
range  of  spatial  frequencies.  Here,  we  restrict  our  analysis  to  the  simple  case  of  a  plane 
wave  incidence,  which  can  be  extended  easily  to  analyze  the  general  case.  Figure  2 
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illustrates  the  progression  of  the  light  through  a  thick  hologram.  This  process  is 
represented  mathematically  by  a  series  of  Mueller  matrices''’^^  that  describe  the 
transformation  of  the  input  Stokes  vector  Si.  The  Stokes  vector  Sfs  of  the  transmitted 
beam  through  the  front  surface  of  the  hologram  is 
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where  Mfs  is  the  Mueller  matrix  for  the  transmission  from  the  front  surface  and  1 1 1  and  ti 
are  the  Fresnel  transmission  coefficients  corresponding  to  the  components  of  linearly 
polarized  light  parallel  and  perpendicular  to  the  plane  of  incidence.  One  can  describe  the 
Mueller  matrix  for  the  exit  surface  of  the  hologram  in  the  same  manner. 

In  order  to  construct  the  Mueller  matrix  for  the  grating  itself,  Mh  in  Figure  2,  we 
describe  the  amplitude  of  the  diffracted  beam  for  parallel  and  perpendicular  incidence  in 
a  manner  analogous  to  the  Fresnel  reflection/transmission  case.  For  our  purposes  here, 
we  need  to  describe  only  the  diffracted  beam  and  not  the  transmitted  component.  The 
relevant  parameters  are  illustrated  in  Figure  3.  U,-  is  the  input  beam  andU^  is  the  output 
beam,  where  and  are  the  respective  polarization  vectors  normal  to  the  direction  of 

propagation.  K  is  the  grating  vector,  which  makes  an  angle  ^  with  the  z  axis.  We  recall 
that  the  amplitude  of  the  diffracted  beam  is  given  by'^: 
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Here,  n'  is  the  index  modulation  depth.  For  the  special  cases  of  polarizations  parallel  and 
perpendicular  to  the  plane  of  incidence,  the  dot  product  in  (2)  can  be  simplified  as: 
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The  Mueller  matrix  for  the  hologram  (Mh)  and  the  transformation  of  the  Stokes  vector  Si 
are  given  by 
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For  the  architecture  proposed  above,  we  must  describe  the  diffraction  amplitudes  from 
two  multiplexed  gratings.  These  gratings  must  be  specially  designed  so  that  they  have 
the  same  Bragg  angle  ^5.  Following  our  previous  analysis,’'*  we  can  design  two 

orthogonal  gratings  such  that  they  each  share  the  same  Bragg  angle.  The  amplitude  of 
the  diffracted  beam  from  the  j”’  grating  is  given  as: 


U;  J  cos(^„y)lo 


sin  (1^0  d)  for 7=1 ,2  where 


(5) 


We  can  evaluate  the  parallel  and  perpendicular  polarization  cases  as  in  (3)  and  obtain  the 
proper  coefficients  U||anduj^. 

The  proposed  architecture  can  now  be  completely  analyzed  using  the  Mueller 
matrices  found  above.  The  Mueller  matrix  for  one  grating  of  the  hologram  is: 


Mt  -  Mgs  •  Mjj .  Mps 


(6) 


where  Mps  is  the  matrix  for  the  front  surface  and  Mgs  is  the  matrix  for  the  exit  surface. 


Using  (1)  and  (4)  in  (6)  we  find  that  the  transformation  of  the  input  Stokes  vector  is  given 
by: 
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Equation  (7)  shows  the  Stokes  vector  that  is  diffracted  from  one  grating.  The  second 
grating  will  produce  an  equation  of  the  same  form  but  with  different  grating  coefficients 
M.  Both  equations  have  the  same  Stokes  vector  as  input,  but  because  the  grating 
parameters  u  differ,  they  will  produce  different  output.  Note  that  using  intensity 
detectors  and  given  the  form  of  the  Mueller  matrix  in  (7),  we  can  only  determine  the  first 
two  Stokes  parameters.  In  order  to  determine  the  other  Stokes  parameters,  one  must  have 
a  Mueller  matrix  with  non-zero  off-diagonal  elements  in  the  third  and  fourth  column. 
This  can  be  achieved  for  example  by  a  rotation  of  the  holographic  grating  about  the  z- 
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axis  and  by  a  similar  rotation  of  the  polarization  axis  of  the  incident  light.  This  simple 
rotation  is  described  by  the  following  Mueller  matrix: 


*  W 


M^ir) 
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0  0  0  1 


(8) 


This  is  the  optimal  ehoiee  because  it  can  be  realized  by  simply  rotating  the  device  when 
in  use.  This  rotation  allows  us  to  find  two  equations  for  the  two  beams  diffracted  by  the 
first  set  of  multiplexed  gratings: 


7,1  =7;  (4  +4  y{QiCOs(2r0-Uisin(2ri))[Ai  j 

^2  =4(^2  +52)  +  (ficos(2r,)-f^/Sin(2yi))(^2  -^2) 


where  the  subscript  1  and  2  denote  the  two  separate  gratings.  The  angle  yi  denotes  the 

angle  of  rotation  for  gratings  1  and  2.  The  next  set  of  equations  comes  from  the  beams 
that  pass  through  the  quarter-wave  plate  and  diffract  from  the  third  and  fourth  gratings. 
The  quarter-wave  plate  adds  a  ;7’/ 2  phase  shift  that  interehanges  the  U  and  V  parameters 
of  the  Stokes  veetor:  cos(£')  -»  eos(£’  +  ;r/2)  =  -sin(£-)  andsin(f)  sin(s  +  7rf2) 

=  cos{s) .  The  second  set  of  equations  is  rotated  an  angle  ^^2  is  thus: 


^13  =  7,-  (4  +  53  ) + {Qi  cos(2;'2 )  +  sin(2r2 ))  (^  “  ^3 ) 
7,4  =  7,-  (^4  +  ^4 )  +  {Qi  cos(2r2)  +  Vi  sin(2x2))  (^4  "  ^4 ) 


(10) 


Because  the  quarter-wave  plate  has  changed  U  into  V,  we  now  have  four 
equations  involving  all  four  Stokes  parameters.  We  are  required  to  measure  the  four 
diffracted  intensities  in  order  to  determine  the  Stokes  parameters.  The  grating  coefficients 
u  and  the  Fresnel  reflection  coefficients  can  be  measured  to  find  the  coefficient  A’s  and 
B’s  which  in  turn  can  be  used  to  pre-weight  the  detector  arrays  for  the  imaging  process. 

To  determine  the  constraints  on  the  device  parameters  that  are  required  for  fully 
determining  the  complete  Stokes  vector,  we  write  equations  (9)  and  (10)  together  in 
matrix  form  as  7,  =  M'4 .  M'  is  the  so-called  measurement  matrix,  Si  is  the  input  Stokes 


vector  to  be  solved  for,  and  It  are  the  four  measured  intensities.  For  our  material 
parameters,  an  example  measurement  matrix  is  given: 

'0.24320  -0.04187  0.23746  0 

0.68596  -0.03614  0.00637  0 

0.28996  -0.04897  0  -0.27773 

0.21291  0.11626  0  0.02050 


M'  = 


(11) 


We  have  verified  numerically  that  this  matrix  determines  the  Stokes  vector  robustly  for  a 
range  of  input  values.  No  attempt  was  made  to  further  optimize  the  condition’^’’^  of  the 
matrix,  which  requires  an  exhaustive  search  through  the  parameter  space,  and  is  subject 
to  further  investigation. 
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To  summarized,  we  have  developed  a  polarization  imaging  architecture  using 
thick  multiplexed  holograms  that  has  many  advantages  over  current  polarimetric  imaging 
techniques.  The  analysis  showing  the  principle  of  operation  explicitly  takes  into  account 
the  polarization  state  of  the  incident  and  diffracted  wave  fields.  Transformation  of  initial 
Stokes  parameters  by  grating  diffraction  is  formulated  by  a  Mueller  matrix  defined  in 
terms  of  diffracted  amplitudes  of  planar  and  perpendicular  polarization  components.  A 
procedure  is  outlined  using  two  sets  of  rotated  orthogonal  gratings  and  a  quarter-wave 
plate  to  compute  all  four  unknown  Stokes  parameters  required  for  polarimetric  imaging. 
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Figure  2:  Interaction  with  hologram  viewed  in  three  sections 
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Figure  3:  Diffraction  from  a  thick,  slanted  grating 
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3.6  Publicity  Articles  Regarding  Super-Parallel  Architecture  and  Slow-Stopped 
Light. 


[SEE  ATTACHMENTS  5,6, 7,8  FOR  SAMPLES] 

4.  TECHNOLOGY  TRANSFER: 

None. 
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